Report for Model 3

— for both jointly and separately

2024-08-19

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.01
## LDevsum       1.02       1.07
## dh0           1.04       1.15
## dh1           1.16       1.55
## dh2           1.14       1.50
## dh3           1.01       1.01
## dl0           1.03       1.13
## dl1           1.05       1.19
## 
## Multivariate psrf
## 
## 1.16

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1211.32188 22616.7719
DIC3 1153.83458 22718.4531
PWAIC 42.23985 259.0377
WAIC 1180.55174 22741.2220

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.05
## LDevsum       1.00       1.00
## dh0           1.04       1.15
## dh1           2.70      10.40
## dh2           1.51       2.48
## dh3           2.93      11.81
## dl0           1.09       1.34
## dl1           1.09       1.35
## 
## Multivariate psrf
## 
## 2.59

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1388.12865 22744.0951
DIC3 1250.68953 22833.8895
PWAIC 91.90702 329.4524
WAIC 1320.22046 22870.8368